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ABSTRACT 

We report on a complete sample of 7 luminous early-type galaxies in the Hubble Ultra Deep Field (UDF) with 
spectroscopic redshifts between 1.39 and 2.47 and to Kab < 23. Using the BzK selection criterion we have 
pre-selected a set of objects over the UDF which fulfill the photometric conditions for being passively evolving 
galaxies at z > 1.4. Low-resolution spectra of these objects have been extracted from the HSTh-ACS grism 
data taken over the UDF by the GRAPES project. Redshift for the 7 galaxies have been identified based on 
the UV feature at rest frame 2640 < A < 2850 A. This feature is mainly due to a combination of Fell, Mgl 
and Mgll absorptions which are characteristic of stellar populations dominated by stars older than ~ 0.5 Gyr. 

The redshift identification and the passively evolving nature of these galaxies is further supported by the photo¬ 
metric redshifts and by the overall spectral energy distribution (SED), with the ultradeep HSTh-ACS/NICMOS 
imaging revealing compact morphologies typical of elliptical/early-type galaxies. From the SED we derive 
stellar masses of ^ IO^^'M© and ages of ^ 1 Gyr. Their space density at <z>= 1.7 appears to be roughly a 
factor of 2-3 smaller than that of their local counterparts, further supporting the notion that such massive and 
old galaxies are already ubiquitous at early cosmic times. Much smaller effective radii are derived for some of 
the objects compared to local massive ellipticals, which may be due to morphological K corrections, evolution, 
or the presence of a central point-like source. Nuclear activity is indeed present in a subset of the galaxies, as 
revealed by them being hard X-ray sources, hinting to AGN activity having played a role in discontinuing star 
formation. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: early-types — galaxies: high- 
redshift — cosmology: observations 


1. INTRODUCTION 

In the local universe, as accurately measured by the SDSS 
(e.g., Baldry et al. 2004), passive early type galaxies with stel¬ 
lar masses larger than 10*dominate the counts of most 
massive galaxies, being a factor of 3 more numerous than late 
type galaxies above this mass threshold. About 1/3 of all the 
stars in galaxies in the local universe are hosted by such ob¬ 
jects (Baldry et al. 2004). The process by which these galax¬ 
ies formed is still unclear. At least some of them may have 
formed at relatively high redshifts, in a process that, for its 
rapidity, is reminiscent of the monolithic collapse scenario 
(Eggen, Lynden-Bell & Sandage 1962). On the other hand, 
others may have been assembled at relatively recent epochs, 
through merging of smaller subunits (Toomre 1977). The for¬ 
mation of massive spheroids is a central problem of current 
theories of galaxy formation (e.g., Loeb & Peebles 2003; Gao 
et al. 2004). 
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It is now well established that up to z ~ 1 a significant 
population of red passively evolving early-type galaxies can 
be found in the field among extremely red objects (EROs; 
Cimatti et al. 2002a; Yan et al. 2004; see McCarthy et al. 
2004 for a review), together with dust reddened systems. The 
EROs are highly clustered (Daddi et al. 2000a; McCarthy et 
al. 2001; Roche et al. 2003; Miyazaki et al. 2003; Brown 
et al. 2005; Georgakakis et al. 2005) as expected for the 
progenitors of local massive ellipticals (Daddi et al. 2001; 
Moustakas & Somerville 2002). The space density of z ~ 1 
early-type galaxies is at least within a factor of 2 of the local 
value (Daddi et al. 2000b; Pozzetti et al. 2003; Bell et al. 
2004; Caputi et al. 2004) implying a modest evolution from 
z = 1, especially for the most massive ones. Their stars ap¬ 
pear fairly old ( ^ 3 Gyr), suggesting even higher formation 
redshifts z ^ 2 (Cimatti et al. 2002a; Treu et al. 2005). 

On the other hand, physically motivated galaxy formation 
models based on hierarchical clustering within the ACDM 
framework (Cole et al. 2001; Kauffmann et al. 1999; 
Somerville et al. 2001) have been so far unable to account 
for the large space density of z > 1 red galaxies (Daddi et al. 
2000b; Smith et al. 2001; Firth et al. 2002; Somerville et 
al. 2004; Glazebrook et al. 2004). These models result in 
widespread merging and associated star-formation activity of 
massive galaxies at relatively low-redshifts, and it is currently 
not fully clear what physical mechanism needs to be involved 
in order to terminate star-formation and produce the red col¬ 
ors of passive galaxies. Feedback processes, e.g., from the 
onset of AGN activity, seem to be a promising tool to achieve 
that (Granato et al. 2001; 2004; Springel et al. 2004). This 
strengthens the necessity to simultaneously trace the forma¬ 
tion of galaxies and AGN activity to understand the link be- 
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tween assembly of stellar mass in galaxies and the growth of 
supermassive black holes (Magorrian et al. 1998; Ferrarese & 
Merritt 2000). 

As little evolution in massive (> 10*'M©) held early-type 
galaxies is detected up to z = 1, it is necessary to push the 
investigation to the highest possible redshifts in order to fur¬ 
ther constrain the formation of early type galaxies. Crucial 
questions to ask are; up to which redshift do passively evolv¬ 
ing early-type galaxies exist ?; which is their space density 
as a function of z ?; what is their clustering ? and what is 
the environment they live in (i.e. cluster vs held) ? Study¬ 
ing the highest redshift massive and passive objects, i.e. those 
presumably closest to the formation or assembly epoch, could 
in principle also reveal useful information to understand the 
physical processes by which these galaxies were formed. 

For almost a decade the unique example of spectroscop¬ 
ically conhrmed high-z passive galaxy beyond z ~ 1.5 has 
been 53W091 at z = 1.55 (Dunlop et al. 1996; Spinrad et 
al. 1997). This objects was preselected for being a radio¬ 
galaxy, hence it was virtually selected over the whole sky. 
Its very red colors suggested fairly old stellar populations, al¬ 
though with some controversies about detailed age (see e.g., 
Nolan et al. 2001 and references therein). The North and 
South Hubble Deep Fields with their extremely deep imag¬ 
ing data prompted searches for z > 1.5 passive galaxy candi¬ 
dates (Treu et al. 1999; Stiavelli et al. 1999; Broadhurst & 
Bowuens 2000; Benitez et al. 1999), however none of such 
objects was eventually spectroscopically conhrmed. That is 
due to the faintness of the candidates in the optical domains 
and because optical spectroscopy from the ground is ham¬ 
pered by the presence of strong OH sky emission lines at the 
wavelengths where the main spectral features are redshifted. 
Now, over helds of a few tens to a few hundreds arcmin^, 
three groups have recently presented the discovery of several 
high-z passive galaxies. Cimatti et al. (2004) reports 4 pas¬ 
sive galaxies spectroscopically conhrmed at 1.6 < z < 1.9 in 
the K20 survey, from a region within the GOODS-South held. 
The ACS imaging revealed their compact early-type galaxy 
like morphologies. McCarthy et al. (2004) have redshifts for 
20 passive objects in 1.3 < z < 2.15 (of which 8 at z > 1.5) 
from the Gemini Deep Deep Survey. Saracco et al. 2005 con- 
hrm 7 very massive old galaxies at 1.3 < z < 1.7 from low- 
resolution near-IR spectroscopy. The latter two studies both 
lack HST morphology information. 

While it is now established that a signihcant population of 
passive galaxies exist up to z ~ 2 at least, the availability of the 
Hubble Ultra Deep Field (UDF) dataset, with ultimately deep 
and high spatial resolution optical and near-IR imaging offers 
a unique possibility for studying these high-z passive galax¬ 
ies in some detail. As a part of our Grism ACS Program for 
Extragalactic Science (GRAPES), we have collected about 50 
HST orbits of ACSH-G800F grism spectroscopy on the UDF 
(Pirzkal et al. 2004, P04 hereafter). As discussed in more de¬ 
tail further in the paper, the HSTh-ACS spectroscopy allows 
us to obtain much higher S/N ratios on the continuum than 
reachable from the ground. In this paper, we have taken ad¬ 
vantage of the GRAPES spectra to present spectroscopic con¬ 
firmation for a sample of 7 early-type galaxies at 1.4 < z < 2.5 
in the UDF, and used the available multiwavelength datasets 
to study their properties. 

The paper is organized as follows: in Sect. |2lwe comment 
on the usefulness of low resolution spectra of passive galaxies 
and a new spectral index is defined for their characterization. 
Sect.|3discusses the color selection of z > 1.4 passive galaxy 


candidates. In Sect. 0] redshift identihcations are presented 
and discussed in detail. Morohological parameters of the 7 
objects are derived in Sect. 0 Stellar population properties, 
including stellar masses and ages, are estimated in Sect. 
The issue of space density evolution of passive early-type 
galaxies to z = 2.5 is discussed in Sect. 0 while in Sect. 0 
the evolution of morphology is discussed. Sect.0presents the 
X-ray detection of two of the galaxies in our sample. Sum¬ 
mary and conclusions are given in Sect. cni 
We assume a Salpeter IMF from 0.1 and 100 M©, and a 
WMAP flat cosmology with HajUm = 0.73,0.27, and h = 
//o[km s-' Mpc-']/100 = 0.71 (Spergelet al. 2003). 


2 . UV SPECTRA OF PASSIVE GALAXIES: THE Mgu,, FEATURE 

The low resolution of the HST grism spectra does not al¬ 
low us to detect individual absorption or emission lines (if 
not for large equivalent widths, Xu et al. 2005) whose iden¬ 
tification generally yields redshift measurements in higher 
resolution spectroscopy. Nevertheless the GRAPES spectra 
are well suited for the identification of broad, low resolu¬ 
tion spectral features like the strong breaks that are com¬ 
monly found among old stellar populations. Fig. 0 shows 
rest frame UV spectra (the region generally accessible to 
GRAPES HSTh-ACS for 1.4 < z < 2.5 objects) from Bruzual 
& Chariot (2003) spectral synthesis models and from Kurucz 
models of stars (Kurucz et al. 1979), smoothed to the typical 
resolution of our data. While the 4000A break and the 3600A 
Balmer break rapidly disappear beyond 1 pm for 1.4, there 

are other strong age dependent features in the UV that are de¬ 
tectable in relatively low resolution spectra. The most promi¬ 
nent one is a bump in the region 2640-2850A that is due to 
the combination of several strong absorption lines, including 
MgII2800 that is the strongest one. The typical shape of this 
region was first used by Spinrad et al. (1996) to measure the 
redshift of the z = 1.55 passive galaxy, and more recently also 
by Cimatti et al. (2004) and McCarthy et al. (2004). We have 
dubbed this feature the Mguv feature, and defined the Mguv 
index as: 


Mguv = 


2xjiiif,dX 


I2525 + 12725 


2825 


( 1 ) 


where the integration ranges (see Fig. 0 are defined in A. 
The Mguv feature is found to be almost independent of the 
spectral resolution for R ^ 50, the typical rest-frame resolu¬ 
tion of the spectra discussed in this paper, and on dust ex¬ 
tinction. Fig. 0 shows the age dependence of the Mguv fea¬ 
ture for passive as well as actively star-forming galaxies. The 
Mguv feature is a key hngerprint of the presence of passively 
evolving stellar populations because it is not present in young 
dust reddened star-forming galaxies (Fig. 0and0, which can 
still produce similarly red overall colors for high dust extinc¬ 
tion. Its secure detection therefore allows one to establish the 
nature of the sources by breaking the age/dusty degeneracy 
among red sources, as well as the measurement of the galaxy 
redshift. 


3 . PASSIVE z> 1.4 GALAXIES: SAMPLE SELECTION 

In order to search for high-redshift, passively evolving 
galaxies among the 10000 galaxies detected in the UDF 
we relied on their peculiar color properties. Candidates of 
z 1.4 passively evolving galaxies were pre-selected follow¬ 
ing the criteria outlined by Daddi et al. (2004b). These were 
calibrated on the complete K20 survey spectra database that 
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Fig. 1.— The mid-UV spectral energy distribution of galaxies from the Bruzual & Chariot (2003) library (left) and of stars from the Kurucz database (right), 
smoothed to ~ 50 A rest frame resolution, similar to the HST spectra of the objects analysed in this paper. Old simple stellar population (SSP) galaxies or low 
mass stars have very red spectra and a characteristic signature with a peak and a dip in the region 2640-2850A, mainly due to a combination of Fell, Mgl and 
Mgll absolutions lines. This Mguv feature appears after a few 100 Myr of passive evolution (visible in A-type stars or later) and gets stronger as age increase 
(notice the 2 order of magnitude UV continuum dimming for SSP galaxies from 0.2 to 2.0 Gyr). Star-forming galaxies with large reddening (dot-dashed line in 
the left panel) can reach similarly red continua in case of very high reddening, but are distinguished from the former because of the lack of any strong spectral 
feature in the UV, like OB-type stars. The vertical lines (both panels) show the 3 windows used for the measurement of the Mguy index (Fq.ITI. 



log Age yr 

Fig. 2.— The age dependence of the Mgyy index (Eq. Q. The 
three solid lines correspond to SSP models for three different metallicities 
(O.4Zo,Z0,2 . 5 Z 0 ; higher metallicity produce a larger Mguv index at fixed 
age for most of the age range). The dotted line is for a constant star-formation 
model with it (B — V) = 1.2. The Mguv index is almost independent on red¬ 
dening and on the resolution of the spectra for R ^ 50, the typical rest-frame 
resolution of our GRAPES spectra. 


includes the z > 1.4 passive and star-forming galaxies de¬ 
scribed in Cimatti et al. (2004) and Daddi et al. (2004a). 
Defining the color difference: BzK = (z — K)ab — {B — z)ab, 


the candidate z 2 passive galaxies can be located with 
BzK < —0.2 n {z — K)ab > 2.5 (Daddi et al. 2004b). 

Objects satisfying the above condition were retained as a 
primary sample. However, based on Fig. 8 (bottom-left dia¬ 
gram) of Daddi et al. (2004b), also galaxies having z — K > 
2.5 and BzK > —0.2 were considered for detailed analysis, 
as young "proto-ellipticals" could be in principle located at 
BzK > —0.2. Moreover, despite the extreme depth of Z?-band 
UDF imaging the errors on B — z colors are found to be large 
for many red objects with z — K > 2.5. 

This pre-selection criteria require B-, z-, and ^f-band imag¬ 
ing. We have used the deep VLTh-ISAAC /Ty-band images^ 
over the GOODS field (Fig.0 (Vandame et al. 2005 in prepa¬ 
ration), and the ultra-deep ACS images^*’ of the UDF for 
the B— and z—bands (Beckwith et al. 2005 in preparation). 
Galaxies were selected over the area of about 12 arcmin^ 
where deep ACS images are available, requiring total mag¬ 
nitudes Kab < 23 as measured by SExtractor MAG_AUTO 
(Bertin & Arnouts 1996). This is ~ 1 mag fainter than reached 
by the K20 survey, while the area is a factor 3 smaller than 
the K20/GOODS area that is contiguous and in part overlap¬ 
ping with the UDF area analysed here (Fig.|3l. 

Fig. 13 show the resulting BzK diagram. In order to in¬ 
crease the accuracy of the color measurements we have not at¬ 
tempted to reduce the ACS resolution to the much poorer 0.5" 
seeing of the K-band data. For the ACS bands we used SEx¬ 
tractor MAG_AUTO measured in consistent apertures defined 
in the z-band (double image mode). SExtractor MAG_AUTO 
measurements were used also for the K-band (single image 
mode). This results in a different physical aperture of the 
magnitudes in the optical versus near-IR bands. However, we 

® Available at "http://www.eso.org/science/goods/releases/20040430/" 
Available at "http://archive.stsci.edu/prepds/udf/udf_hlsp.html" 
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Fig. 3.— The relative layout of the UDF, K20 and GOODS-South regions. 
The background image is the GOODS K-band mosaic used in this paper. 


verified that the resulting z — K was not biased (at least within 
^ 0.05 mags on average) with respect to seeing matched z — ZT 
magnitudes for objects in common to the catalog of Daddi et 
al. (2004b). There are 20 objects with z — K > 2.5 of which 7 
have also BzK < —0.2. 

4. ANALYSIS: REDSHIFTS AND SPECTRA 

To measure the redshifts of the selected candidates, we have 
used all the available information including high quality HST 
low resolution optical spectra and the multicolor photomet¬ 
ric data available for the UDF, as discussed in the next two 
sections. 

4.1. ACS-HST low-resolution spectroscopy 

The ACS-HST grism 800L data, obtained as part of the 
GRAPES project, cover the wavelength range 5000-1 lOOOA 
with maximum efficiency around 7000-8000A. The spectral 
dispersion is about 40A/pixel. The low sky background at 
these wavelengths from space allow HST grism data to be 
taken in slitless mode, so that all UDF objects are simulta¬ 
neously observed in the GRAPES data. The effective spectral 
resolution depends on the spatial size of the objects. For un¬ 
resolved sources the spectral dispersion translates into a res¬ 
olution of R « 100 at 8000A, while extended sources have 
lower resolution spectra. Thanks to the absence of OH sky 
emission lines and the exquisite spatial resolution of HST, the 
GRAPES low-resolution spectroscopy allows us to reach rel¬ 
atively high S/N ratios on the continuum. For compact galax¬ 
ies, like the ones considered in this paper, one gets S/N 10 
at 8000A over a lOOA region for Iab ~ 25.5, for the lO^s 
integration time of the GRAPES data. This is significantly 
better than what can be reached from the ground with red- 
optimized CCDs on 8-10m class telescopes with similar inte¬ 
gration times, except for the narrow spectral windows free of 
OH sky emission lines. 



Fig. 4.— The BzK diagram for galaxies with Kab < 23 in the UDF. Objects 
with z — K > 2.5 were retained as possible high-z old galaxies candidates. 
The seven passive early-type galaxies with proposed redshift 1.4 $ z $ 2.5 
are shown with large symbols. Five occupy the region with BzK < —0.2 as 
expected (Daddi et al. 2004b), while two bluer objects are also identified. 


Full details of GRAPES grism data reduction and calibra¬ 
tion are given in P04. As slitless spectroscopy can imply spec¬ 
tra superposition from neighboring sources, the data were ob¬ 
tained at five independent position angles to minimize this ef¬ 
fect. Narrow extraction windows (see P04) were used to max¬ 
imize the S/N ratio of the resulting spectra. The 5 epoch spec¬ 
tra were co-added, avoiding regions contaminated by neigh¬ 
boring objects, as described in P04. 

4.2. Spectral analysis 

The 20 red galaxies selected as described in Sect. 0 were 
inspected by looking for the Mguv feature described above, 
and to their overall spectral shape. Some of the spectra had too 
low S/N ratio to be useful because the objects are too faint in 
the optical. Seven objects were retained as likely z ^ 1.4 pas¬ 
sive objects identifications, five of these having BzK < —0.2. 
Spectra for these seven galaxies are shown in Fig.|3 The five 
single epoch spectra of each object were inspected to verify 
that the detected features were persistent and not due to obvi¬ 
ous spurious artifacts. 

In order to objectively identify the redshift of the ob¬ 
jects, and to exclude possible features misidentification, we 
cross correlated the spectra with galaxy templates from the 
Bruzual & Chariot (2003) library before definitively assign¬ 
ing redshifts. The model templates were convolved with the 
line spread functions computed from the light profiles of the 
galaxies to reproduce the actual spectral resolution of the data. 
Note that this fitting approach is rather conservative because 
most of the continuum region of the spectra with no detected 
features have no strong effect on the fitting and reduce the im¬ 
pact of the regions where the actual features are, i.e. it dilutes 
the signal by giving equal weight to all data in the spectral 
range. 

Two classes of models were adopted for the spectral fits in 
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Table 1. Photometric properties, SED and spectral fitting results 


ID 

RA 

DEC 

z"" 

AB 

Kfoi 

AB 

Photometry 
B^z Z-K 
AB AB 

R-K 

Vega 

J-K 

Vega 

SED fitting 

2 

Zphot Xr.old 

xli 

Spectral fitting 
Zspec Xr,old 

xli 

8238 

03:32:36.9 

-27:46:28.5 

23.80 

21.22 

4.48 

2.58 

6.43 

2.01 

1.26 ±0.04 

1.1 

2.8 

1.39±0.01 

0.59 

0.61 

4950 

03:32:30.0 

-27:47:26.8 

22.56 

19.90 

2.12 

2.64 

5.58 

1.96 

1.57 ±0.03 

0.4 

7.2 

1.55±0.01 

0.83 

0.79 

1025 

03:32:43.0 

-27:48:45.1 

24.20 

21.48 

1.83 

2.72 

5.41 

1.87 

1.72 ±0.04 

0.6 

5.1 

1.73±0.01 

0.40 

0.54 

3523 

03:32:33.7 

-27:47:51.1 

24.34 

21.68 

3.41 

2.66 

5.84 

1.79 

1.72±0.04 

2.2 

9.0 

1.76 ±0.02 

0.70“ 

0.64 

3650 

03:32:38.1 

-27:47:49.8 

24.18 

21.36 

3.94 

2.82 

6.00 

1.84 

1.90 ±0.03 

2.2 

9.7 

1.91±0.01 

0.76 

0.98 

3574 

03:32:39.1 

-27:47:51.6 

25.29 

22.39 

3.24 

2.90 

6.29 

1.89 

1.94 ±0.04 

3.0 

11.3 

1.98 ±0.02 

0.41 

0.46 

1446 

03:32:39.2 

-27:48:32.4 

26.06 

22.80 

>3.80 

3.26 

5.42 

3.32 

2.85 ±0.05 

2.1 

5.8 

2.47 ± 0.02 

0.58 

0.71 


N ote. — IDs for UDF galaxies are those from the publicly available catalog release VI, available at "ftp://udf.eso.org/archive/pub/udf/acs-wfc/h_udf_wfc_Vl_cat.txt". Note that the total magnitudes were derived from Sdrsic profile fitting, as described in 
Sect^^and are brighter than the Kron total magnitude measured by SExtractorby ~ 0.4 mags, on average. Conversion factors from Vega to AB magnitudes are (-0.09,0.22, 0.94,1.87) for (B. R, J and K), respectively. Photometric redshift errors are at the 
68% level as computed by hyperz- Spectro.scopic redshift errors are the formal 68% range on the fit determined following Avni (1976) for the case of one interesting parameter (redshift), added in quadrature to systematic and random errors in the wavelength 
calibration that are about 20A (Pasquali et al. 2003: P04) and produce an additional uncertainty of o-j % 0.01. The spectral fit for object #3523 has a slightly lower minimum for z 1.1, a redshift that is however fully inconsistent with the photometric 
SED. 


order to bracket the general cases of old and passive galax¬ 
ies versus strongly dust reddened galaxies. For old galaxies 
we considered simple stellar populations (SSP) and a model 
with exponentially declining SFR with timescale r = 0.3 Gyr. 
Ages less than the age of the universe at each redshift were 
required and a maximum reddening of E{B — V) < 0.2 was 
allowed. We considered solar and 40% solar metallicity. For 
dusty galaxies we used constant star-formation rate models 
with unlimited*^ dust reddening and solar metallicity. 

In all seven cases, except for #3523, the models for old pas¬ 
sive galaxies converge to the solution guessed by eye, with 
the Mguv feature correctly identified. In most cases dusty 
models provide a worse fit to the data, however, differences 
are not large and all reduced are below 1 and acceptable. 
The minima of reduced being all below 1 might actually 
suggest that some significant degree of correlation (e.g., due 
to the reduction process, P04) is present in the data. Results 
of the spectral fitting are summarized in Tabled where in all 
cases the Zspec is quoted for the best fitting old/passive galaxy 
model, as justified from the results of the next section. De¬ 
tailed description of redsh ift assignment for individual objects 
is described in Sect. EH following the analysis of the SEDs 
and photometric redshifts. 

4.3. Spectral energy distributions and photometric redshifts 

To further constrain the redshift identifications we used the 
extremely high-quality imaging available on the UDF area. 
This includes the BViz ACS imaging and the J and H NIC- 
MOS imaging. ISAAC J and K band data were also used, the 
J-band filter of ISAAC being much narrower than the NIC- 
MOS one, as well as V and R from FORS2 that have cen¬ 
tral wavelengths that are quite different from the ones of the 
ACS bands. To match ground based and HST datasets we pro¬ 
ceeded as discussed in Sect. 0 It was checked that over the 
common range the photometric SEDs and the HST low reso¬ 
lution spectra were consistent. The SED of the seven z > 1.4 
early-type galaxies candidates are shown if Fig.|6l The same 
range of models parameters as for the spectral fitting was 
used, corresponding to the two classes of old/passive versus 
dusty models. The hyperz code (Bolzonella et al. 2000) was 
used for the model fitting and to derive photometric redshifts. 
A lower limit of 0.05 mags was used for the photometric er¬ 
rors in all the bands, in order to account for photometric ze- 
ropoint uncertainties and possible residual uncertainties in the 

** practically, for the intrinsically bluest models, E{B — V) iS 1.5 is more 
than sufficient to provide SEDs as red or redder than the data at any redshifts 


match of the different datasets. 

Results of the SED fitting and photometric measurements 
are summarized in Tabled SED fitting in all cases result in 
significantly better fits for models of old/passive galaxies with 
respect to those for dusty objects, rejecting dusty models at 
more than 3 (t confidence levels (Avni 1976) in all cases. Best 
fits have generally acceptable reduced and indeed look rea¬ 
sonable by eye inspection. The agreement between the photo¬ 
metric redshifts and the spectroscopic redshifts, derived from 
the spectral fitting with old stellar populations, is very good 
with 6 out of 7 objects having \zspec — Zphot I ^0.1. 

4.4. Redshifts and redshifts quality classes for individual 

objects 

Based on the results of SED and spectral fitting, we as¬ 
signed redshifts and redshifts quality classes (A for good- 
confidence redshifts, B for less secure ones) for individual 
objects. A quality "A" redshift is assigned when the Mguv 
feature appears well detected, and when the photometric and 
spectroscopic redshifts are in good agreement (i.e., within 
Az = 0.1). 

• #8238: z= 1.39, Class = B. 

The SED and spectral fitting agree within Az = 0.13, the spec¬ 
trum is consistent with the presence of the 4000A break and 
of the Mguv feature. 

• #4950: z= 1.55, Class = A. 

The SED and spectral fitting agree very closely with Az = 
0.02, the spectrum clearly shows the Mguv feature. The rise 
toward the UV below 6000A is consistently present also in 
the photometry. This object is also part of the K20 survey 
(Cimatti et al. 2002b) with a measured redshift of z = 1.553. 
A weak [OII]A3727 emission line is detected in the K20 sur¬ 
vey spectrum. 

• #1025: z= 1.73, Class = A. 

The SED and spectral fitting agree extremely closely with 
Az = 0.01. The spectrum, that is fairly blue especially at 
A < 2500A, shows the characteristic Mguv feature. The up¬ 
turn in the blue spectra is consistent with what expected for 
A-type stars (Fig. Q. This source is detected in the X-r^s 
(Giacconi et al. 2002: Alexander et al. 2003; see Sect. |^. 
Zheng et al. (2004) propose 1.47 < Zphot < 1-58, reasonably 
consistent with our spectroscopic redshift. 

• #3523: z= 1.76, Class = B. 

The spectrum would suggest z ~ 1.1 as the best fitting old 
galaxy solution by placing the 4000A break around 9000A 
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Fig. 5.— The GRAPES HST+ACS spectra of the seven passively evolving galaxies with 1.39 < z < 2.47. In each panel the observed spectra are fitted with 
Bruzual & Chariot (2003) models for old galaxies (red) and main sequence stars from the Kuruez library (green). Proposed redshifts and UDF IDs are labeled on 
each panel. Both rest-frame (top x-axis) and observed (bottom x-axis) wavelengths are labeled. Color inserts for each galaxies are 5'^ x 5'^ (except for #1025 that 
has a smaller color image because close to the edge of the UDF field). 


and the Call H&K lines in the dip at 8700A. The above solu¬ 
tion is not consistent with the photometric SED, as at z Ri 1.1 
no reasonably good fit could be obtained. A secondary so¬ 
lution from spectral fitting, with nearly the same of the 
best one, is found at z = 1.76, in very good agreement with 
the photometric redshift within Az = 0.04. We propose the 
latter solution as the spectroscopic redshift for this galaxy. 
The Mgijv feature is not evident, although the spectrum is 
consistent at least with the presence of a 2900A break. The 
proposed redshift thus relies strongly on the overall SED and 


spectral fitting. The absorption like spectral features at 8900A 
and 9200A are to be considered spurious noise features for 
our proposed redshift. Yan et al. (2004) suggest this is an 
old galaxy at Zphot = 1 -6, in good agreement with the redshift 
proposed here. 

• #3650: z= 1.91, Class = A. 

The SED and spectral fitting agree extremely closely with 
Az = 0.01. The Mguy feature is very evident. 

• #3574: z= 1.98, Class = A. 

The SED and spectral fitting agree closely with Az = 0.04. 
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Fig. 6 .— The observed spectral energy distributions of the seven, z ^ 1.4 early-type galaxies. The imaging data shown in the plot are, from blue to red: F435W, 
V, F606W, R, F775W, F850LP, FI lOW, J, F160W, K. The best fit Bruzual & Chariot (2003) SSP models overplotted here are consistent with those used in Fig. 5. 
Redshift is increasing from bottom to top. The upper limits shown are at Itr. 


The spectrum shows the Mguv feature. Yan et al. (2004) pro¬ 
pose this is an old galaxy at Zphot = 1 -9, in excellent agreement 
with our spectroscopic redshift. 

• #1446: z = 2.47, Class = B. 

The SED and spectral fitting are in fair agreement at best with 
Az ~ 0.4. The spectrum has a strong feature at 9200A that 
in the spectral fitting with old models is identified with the 
Mguv feature. Emission line identification algorithms de¬ 
scribed in Xu et al. (2005), however, pick up the feature as 
a possible emission line. If the feature is an emission line we 
find best fitting A ~ 9230A, EWHM~ 160A and EW~ 150A. 
The large EWHM and EW would suggest an AGN emission 
line. This would be possible as the object is an hard X-ray 
source (Giacconi et al. 2002; Alexander et al. 2004), although 
no strong emission line is detected in the EORS spectra de¬ 
scribed by Szokoly et al. (2004). We tried to fit its photomet¬ 
ric SED with reddened QSO/AGN templates but the fit still 
significantly prefers old/passive models. As the SED drops in 
the B-band, the only plausible AGN line identification would 
be [CIII]A1909 at z = 3.83, with the Lyman-break producing 
the B-band drop. Zheng et al. (2004) propose a photometric 
redshift 4.13 < Zphot < 4.35, not far from that. Eor z = 3.83 
one would expect to detect an even stronger CIVA1550 emis¬ 


sion at about 7500A that instead is not present. Eor z = 3.83 
it would be difficult also to explain the overall shape of the 
SED with, e.g., the break between the J- and H-bands. Chen 
& Marzke (2004) propose this object is a dusty galaxy with 
Zphot = 3.43, not consistently with our analysis, while Yan et 
al. (2004) propose this is an old galaxy at Zphot = 2.8. 

The spectral and photometric properties show that we can 
distinguish two classes of objects among our sources. Objects 
#4950 and #1025 have in fact significant B-band excess flux 
with respect to the others (see Eig.l^and Eig.|^ suggesting the 
presence of hotter star s. Th ese two objects will be discussed 
in more details in Sect. 16.21 

4.5. Discussion of redshift identifications 

The proposed redshifts are based on low resolution spectra 
that, we recall, do not allow to detect individual absorption or 
emission lines as ordinarily done for redshift identifications 
of faint galaxy spectroscopy. Therefore, we cannot exclude 
with full confidence that, in a minority of the cases, our iden¬ 
tifications might be somehow mistaken. We believe, how¬ 
ever, that identifications are generally reliable, because: (1) 
the spectral fitting routine generally converge to the solution 
guessed by visual inspection of the spectra, when using tern- 
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Fig. 7.— The Pozzetti & Mannucci (2000) diagnostic diagram with the 
R — K versus J — K colors of our sample galaxies. The dotted lines show the 
R — K > 53 limit of validity of the criterion and the division between old and 
dusty star-forming EROs for z 5 2. 


plates for passive galaxies; (2) the photometric redshift esti¬ 
mated with the high-quality UDF data are in very good agree¬ 
ment with the proposed spectroscopic redshifts with median 
\zspec — Zphot\ = 0.04 only; (3) there is an overall consistent 
picture as signatures of old stellar populations are seen in the 
spectra, the SEDs strongly argue for passive stellar popula¬ 
tions and morphologies of early-type galaxies are recovered, 
as discussed in the next section. 

As an additional test, we verihed the location of the seven 
sources in the Pozzetti & Mannucci (2000) diagram, Eig. 0 
The R — J versus J — K colors of the sources are consis¬ 
tent with those expected for passively evolving galaxies up 
to z ^2. The only due exception is object #1446 with much 
redder J — K color and proposed z = 2.47, a redshift at which 
the above diagnostic may not apply. We verified that using the 
criterion proposed by Pozzetti & Mannucci (2000) for z> 2, 
involving the J-, H- and K- bands, also the object #1446 would 
be classified as a passive red galaxy. 

We also verified, for further overall consistency, that best 
fitting models of the SEDs produce good fits to the HST grism 
spectra. 

4.6. Sample completeness 

The completeness of the z sample of passive galaxies 
was investigated. We checked if additional passive z ^ 1.4 
galaxies could be found among bluer objects with 2.2 < 
z — K < 2.5. It may be expected that toward the lowest 
z ~ 1.4 range of redshifts some early type galaxy could have 
z — K< 2.5 (either intrinsically or because of the errors in the 
photometry). No convincing case could be found. Among 
the remaining unidentified galaxies with z — K >2.5 there are 
3 more objects with compact morphology that are not in our 
sample with proposed spectroscopic redshift. One of these is 
#869, which is a z = 3.064 type 2 QSO (Szokoly et al. 2004), 
with no distinctive feature in our GRAPES spectrum. The 


other two are UDF #8363 and #5056, and have noisy spectra 
with no feature identifiable and photometric redshifts 2.30 and 
1.43, respectively. These might be additional passive z > 1.4 
objects. Their BzK colors are consistent with the BzK < —0.2 
region (object #8363 is the reddest one in z — with a lower 
limit to B —z, Fig.|3. 

There are 10 more galaxies with z — K > 2.5 that remain 
unidentified in GRAPES spectra, and have irregular, diffuse 
and sometimes merging-like morphologies, reminiscent of the 
z 2 star-forming galaxies (having also similar colors and 
magnitudes) shown in Daddi et al. (2004a,b). Only one of 
these has formally BzK < —0.2. The fact that we cannot 
detect features of old stars in these objects could be an ob¬ 
servational bias as the S/N of the spectra for these extended 
and diffuse galaxies is generally low. However, they might 
be more likely genuine dust reddened star-forming galaxies. 
An analysis of the stacked X-ray and radio emission of galax¬ 
ies with similarly z — K colors and having BzK < —0.2 in the 
K20 survey suggests that these objects are mainly vigorous 
dust reddened starbursts at z > 1.4 (Daddi et al. 2004a,b). 

In summary, there are at most only two reasonable addi¬ 
tional candidates passive z > 1.4 galaxies, that are compact 
and have the expected BzK colors. These are less luminous 
than the ones for which we propose a spectroscopic redshift. 
Additional passive objects might still be present among the 
galaxies with irregular morphology, although no object with 
similar properties is currently known both at high and lower 
redshift. 

4.7. Redshift clustering 

We mention in passing the few evidences of redshift clus¬ 
tering that can be drawn from the redshifts of our sample. For 
these objects, the correlation length ro might be comparable 
or even higher than for typical EROs at z ~ 1 (ro ^ 10 
Mpc; Daddi et al. 2001; McCarthy et al. 2001; Brown et 
al. 2005), as these could be the truly first collapsed massive 
galaxies. This should translate in significant pairings between 
ours and to other galaxies, e.g. to those in the K20 region 
(Fig.|3. In fact, the z = 1.91 galaxy is at the same redshift of 
the Cimatti et al. (2004) z = 1.903 elliptical and also to that of 
the z = 1.901 near-IR bright galaxy in Daddi et al. (2004b). 
The z = 1.73 and z = 1.76 galaxies in our sample may actually 
be at the same redshift given the errors, presumably z = 1.73 
where two other near-IR bright massive star-forming galaxies 
also lie (Daddi et al. 2004b). Also at z = 1.55 and z = 1.39 
other massive star-forming galaxies are found in the K20 sur¬ 
vey. Similarly, 3 out of 4 of the z > 1.5 early-type galaxies 
by Cimatti et al. (2004) are at z 1.61, and strong pairing is 
observed also in the McCarthy et al. (2004) sample. 

5 . ANALYSIS: MORPHOLOGY 

We inspected by eye the morphological appearance of the 
7 galaxies with proposed spectroscopic confirmation as old 
high-z stellar populations. Fig. |3 show the available HST 
imaging (see also the color images of the objects in Fig. |3i. 
All objects look very compact and regular. Object #4950 ap¬ 
pears a bulge dominated nearly face-on spiral, i.e. an Sa, with 
perhaps evidence of merging or cannibalism in the blue bands. 
Object #1446 is the most elongated one and might be an early 
spiral or an SO galaxy. The other objects are visually classi¬ 
fied as E/SO systems. In order to provide a more quantitative 
estimate of the morphology we used CAS and Sersic profile 
htting. 
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Fig. 8. — Multiband HST imaging of the high-z early-type galaxies. Each image is 2” on a side. Within the ACS imaging, as within the NICMOS imaging 
but on a different scale, the same range of fluxes {f\) are plotted in a logarithmic scale, to provide direct comparison with Fig. I^and among different objects or 
different bands. The F435W band appears the noisiest in f\, where however only object #1446 is not significantly detected. The NICMOS images allow in some 
cases to see the extended, low-surface brightness wings of the steep de Vauculeurs like profiles. Objects #1025 and #4950 are outside of the region covered by 
the NICMOS imaging. 


5.1. CAS analysis 

The concentration and asymmetry of all the galaxies hav¬ 
ing z—^T >2.5 were measured in the /-band, corresponding to 
about 3000A rest-frame for z^2 (similar results would be ob¬ 
tained in the z-band). Parameters definitions consistent with 
those given in Conselice (2003) were used. Fig.shows the 
resulting measurements. Objects with proposed spectroscopic 


confirmation as z > 1.4 passive galaxies tend to have high con¬ 
centrations and low asymmetry (see also Tab. Ell. In general, 
these objects occupy a range of values consistent with those 
of early-type galaxies; for example all have concentrations 
C > 2.6, a value marking the boundary between early- and 
late-type galaxies in the local universe (Strateva et al. 2001; 
see also Conselice 2003). The asymmetries are generally low: 

0.2 in all cases. 
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Table 2. Morphological parameters 


ID 

Type 

/-band 

C A 

n 

re pix 

Z-band 
fe Lpc 

b/a 

n 

re pix 

/-band 
re kpc 

b/a 

8238 

E/SO 

3.09 

0.11 

8.2±1.5 

11±4.7 

2 .8±1.2 

0.89 

5.8±1.0 

4.7±1.I 

1.2±0.3 

0.88 

4950 

Sa+M? 

3.57 

0.18 

4.3±0.4 

22 ±3.7 

5.6±0.9 

0.60 

- 

- 

- 

- 

1025 

E/SO 

3.59 

0.12 

4.2±0.5 

2.9 ±0.3 

0.74±0.1 

0.85 

5.0±0.7 

2.1 ±0.3 

0.54 ±0.1 

0.85 

3523 

E/SO 

4.16 

0.19 

9.0±1.6 

11±4.3 

2 .8±1.1 

0.74 

9.7 ±2.0 

16±8.0 

4.0 ±2.0 

0.71 

3650 

E/SO 

2.87 

0.12 

4.7±0.6 

3.1±0.3 

0.79 ±0.08 

0.74 

5.4±0.8 

3.2±0.5 

0.81±0.13 

0.72 

3574 

E/SO 

2.72 

0.12 

2.9±0.3 

2.5±0.2 

0.63 ±0.05 

0.45 

2.8±0.4 

2.4 ±0.4 

0.61 ±0.1 

0.39 

1446 

SO/Sa 

2.81 

0.16 

0 .8±0.2 

3.1±0.5 

0.76±0.12 

0.35 

1.4±0.3 

3.9±0.8 

0.96 ±0.20 

0.26 


5.2. Sersic profile fitting 

We used GALFIT (Peng et al. 2002) to model the light 
profile of the 7 objects with proposed redshift confirmation. 
Both /-band and z-band analyses were performed for a cross 
check of the results and to test for possible color trends. The 
PSF was derived for each band by averaging stellar objects in 
the field (Pirzkal et al. 2005). All the z ^ 1.4 passive galaxies 
are resolved. We fitted Sersic (1968) profiles of the form: 

/L(r) =/Lg ^ (2) 

where is the effective radius, gte is the effective surface 
brightness, n is the free Sersic index, and k is determined 
from n in order for half of the inte^ated flux to be within 
Vg. Results are summarized in Table |2] The Sersic profile fit¬ 
ting also allowed us to derive the total magnitudes, as listed 
in Tabled The quoted uncertainties are purely statistical and 
derived by GAFFIT on the basis of the image noise. Notice 
that the errors on the parameters are correlated. For example, 
fitting with a larger n generally implies larger Vg and brighter 
total magnitudes. 

The Sersic index n is generally found consistent within the 
estimated noise between the i- and z-bands. Derived values 
are in the range of 2.8-9 with many of them close to the 
de Vauculeurs value of n = 4, typical of E/SO systems, with 
the exception of the object with the highest estimated redshift 
(#1446) having n sa 1, that is consistent with the presence of 
an exponential disk and may be an early-type spiral or per¬ 
haps an SO galaxy. This galaxy is reminiscent of the z = 2.5 
evolved disk galaxy discussed by Stockton et al. (2004). For 
object #4950 the Sersic fit leaves strong residuals, correspond¬ 
ing to the features visible in the B-band (see Fig. |^. Never¬ 
theless, the z-band light appears to be dominated by a regular 
spheroid with n ~ 4. A fit in the /-band was not attempted for 
this object. 

The effective radii in the /- and z-bands are also in most 
cases consistent, suggesting the lack of strong color gradients 
within the galaxies. Some of the typical effective radii are 
very small, of order rg ^ 2^ ACS-UDF pixels, correspond¬ 
ing to 0.06-0.12" for the 0.03"/pix scale of the drizzled UDF 
images, or to a physical size 1 kpc (proper length). We will 
come back to these issues in Sect.|3 The axial ratios bja are 
also consistent between the two bands, generally quite close 
to unity and consistent with the range expected for E/SO sys¬ 
tems. 

The proposed reds hift identifications result in a reasonably 
complete (Sect. 14.6> . magnitude limited sample of z > 1.4 
galaxies representative of old stars dominated, morpholog¬ 
ically established early type galaxies. These high-redshift 
spheroids appear to be already dynamically relaxed, with the 
possible exception of object #4950. 


0.6 

0.5 

0.4 

< 
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0.2 

0.1 
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Fig. 9. — The concentration versus asymmetry parameters in the /-band for 
the 20 galaxies with z — K > 2.5 and Kab < 23 in the UDF. Large circles are 
for the seven z ^ 1.4 early-type galaxies. A cross around a symbol is for ob¬ 
jects classified as "compact" by visual inspection, as opposed to the remain¬ 
ing objects (shown without crosses) that are classified as "irregular/merger". 
The three compact objects that are not i n ou r list of passively evolving galax¬ 
ies are described in more detail in Sect. l4.6l 


6. RESULTS: STELLAR POPULATIONS OF z> 1.4 EARLY-TYPE 
GALAXIES 

In order to derive the characteristic stellar population prop¬ 
erties of the seven z > 1.4 galaxies we compared the photo¬ 
metric SEDs (with deep imag ing data in 10 bands from B to 
K, as described in Sect. EJ see also Fig. |5J to stellar pop¬ 
ulation synthesis models from the Bruzual & Chariot (2003) 
library, at the (fixed) spectroscopic redshift of each source. 
We base our physical parameter estimates mainly on the pho¬ 
tometric SEDs rather than on the spectra because of the much 
wider wavelength range spanned by the former. In particu¬ 
lar we would expect the rest-frame UV photometry, including 
the very deep B-band images, to strongly constrain the pres¬ 
ence or lack thereof of young stellar populations. The near-IR 

bands provide instead a census of the past star-formation. _ 

A more extended library than the one adopted in Sect. 14.31 
was used in this case, including models with 0.2, 0.4, 1.0 and 
2.5 solar metallicity. Reddening is restricted to be £ (B — V) < 
0.2 (we use a Calzetti et al. 2000 law). Some small amount of 
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Table 3. Stellar populations properties 


ID 

Zspec 

Class 

Xlm 

nxlJ 

% 

Mb 

Vega 

U-Bi-est 

Vega 

B-V,est 

Vega 

M* 

10 '‘Mq 

M^/Lb 

0 

A-gepass 

Gyr 

■^pass 

‘^type 

Mguv 

8238 

1.39 

B 

1.5 

13 

-20.98 

0.23 

0.91 

1.0-2.4 

2.5-6.2 

0.5^.5 

> 1.6 

F5 

1.14±0.19 

4950 

1.55 

A 

1.2 

30 

-22.86 

0.19 

0.76 

2.9-13 

1.3-3.4 

0.6-2.5 

1.8-3.7 

FO 

1.17±0.04 

1025 

1.73 

A 

0.7 

68 

-21.76 

0.19 

0.76 

1 .0-1.8 

1 .3-2.2 

0 .6-1.2 

2 .1-2.6 

A3 

1.14±0.06 

3523 

1.76 

B 

1.7 

8.4 

-21.73 

0.20 

0.67 

1 .0-1.5 

1.3-1.8 

0.5-1.2 

2 .1-2.6 

F2 

0.94 ±0.09 

3650 

1.91 

A 

1.7 

7.8 

-22.32 

0.20 

0.67 

1 .3-2.0 

1.0-1.5 

0 .6-1.6 

2.3-3.4 

FO 

1.25 ±0.05 

3574 

1.98 

A 

2.8 

0.2 

-21.41 

0.16 

0.68 

0.5-0.9 

0.9-1.6 

0.5-1.2 

2.3-2.9 

FO 

1.5±0.2 

1446 

2.47 

B 

10 

0 

-21.62 

0.20 

0.59 

0.7-1.1 

1.0-1.7 

0.3-0.7 

2.8-3.2 

F5 

2 .8±0.6 


Note. — Xrm reduced the best fit. Range.s for A/# and Agepass (and thus forM */Lg and zpass) are given \ 


the 95% confidence level. U —B and B — V colors are derived from the best fitting model. 


dust reddening are sometimes found also among local early- 
type galaxies (e.g., Goudfrooij & de Jong 1996), and could be 
more often present at these high redshifts as we are observ¬ 
ing more closely to the star-formation epoch. For the star- 
formation history we adopted SSP models and exponentially 
decreasing models with t = 0.1, 0.3, 1.0 Gyr. In addition 
we considered models with step-wise star-formation history, 
with constant SFR lasting for 0.1, 0.3, 1 and 2 Gyr, followed 
by a period of passive evolution. Only ages less than the age 
of the Universe at the fitting redshifts were allowed. The pa¬ 
rameters we are mainly interested in are the age and stellar 
mass of the galaxies and we determined 95% ranges for these 
two parameters following Avni (1976) and marginalizing over 
E{B — V), metallicity and SF history (Table^. We also fitted 
main sequence stars templates from the Kuruez library to the 
HST spectra of all galaxies (Fig.|3. The Mguv indexes were 
measured. Indeed, in most cases they are significantly larger 
than 1, as expected for old stellar populations. 

All the results are shown in Tabled Best fits at the spec¬ 
troscopic redshifts are shown for all objects in Fig. The 
reduced Xm generally low, with large associated probabil¬ 
ities, and the best fits to the photometric SEDs (Fig. |6|i look 
generally very good. An exception is object #1446 with large 
Xm at the tentative spectroscopic redshift z = 2.47. Still, its 
SED looks reasonably well reproduced at the eye inspection 
(Fig- EJ. In part the bad Xm Fe due to the presence of 
AGN continuum emission starting to contribute to the near- 
IR bands fluxes. 

The derived masses are generally large, in the range 0.5- 
7.3 x 1O"M0. Notice that these values were derived using 
the total flux of these galaxies estimated from Sersic profile 
fittings. The use of IGon magnitudes derived by SExtrac- 
tor would have implied an underestimation of the masses by 
~ 40% on average. Absolute magnitudes and colors were de¬ 
rived for our targets by using the best fitting models for com¬ 
puting K-corrections. The derived B-band rest-frame stellar 
mass to light ratios are typically about 1-2 in solar units, as 
compared to the dynamically estimated total (thus including 
the contribution of dark matter) mass to light ratios of up to 6 
for local early-type galaxies (van der Marel 1991). Account¬ 
ing for the fact that the adopted Salpeter IME all the way down 
to 0.1 Mq may be over-estimating the stellar masses (and mass 
to light ratios) by at least 30%, and that dust reddening in the 
fit can affect by up to 1 magnitude the B-band rest-frame, this 
is consistent with a dimming of a few magnitudes from pas¬ 
sive evolution to z = 0 as computed from aging the best-fitting 
Bruzual & Chariot models. The inferred stellar mass to light 
ratios at z = 1.4-2 appear also consistent with the trend estab¬ 
lished for cluster early-type galaxies up to z ~ 1.3 on the basis 


of fundamental plane studies (van Dokkum et al. 2004; Ren- 
zini 2004). The 95% ranges on fitted properties are not very 
narrow in most cases. Availability of Spitzer photometry pro¬ 
viding access to the near-IR rest-frame could allow to further 
improve these estimates. However, the uncertainties about the 
importance and modeling of AGB stars contribution at those 
wavelengths (Maraston 2005) could further complicate the is¬ 
sue, our sample being perhaps the ideal one to verify the pre¬ 
dictions of different models. A future publication will address 
these points. 

6.1. Formation and ages of z .^1.4 early-type galaxies 

In order to characterize the evolutionary status of these ob¬ 
jects, we estimated the quantity Agepass as the time elapsed 
since the onset of the passive evolution (i.e. since the end 
of the last major burst of star-formation). This is defined as 
the difference between the model age and the duration of the 
star-formation phase, assumed to be equal to r for exponential 
SE histories, zero for SSP models, and the actual duration of 
the burst for truncated star-formation models. Tabled shows 
that Agepass estimates are typically in the range of 0.5-1.5 
Gyr, suggesting that the onset of passive evolution started at 
Zpass ~ 2-3.5 in most cases. 

It is interesting to investigate what types of SE histories 
yield acceptable fits to the SEDs, as this could constrain the 
SE timescales in early-type galaxies. Limiting to the 5 ob¬ 
jects with no B-band excess flux, we find that the SSP and 
T = 0.1 Gyr models provide much better fits to the SEDs than 
models with r =0.3 or 1 Gyr. Exponentially declining star- 
formation histories with r >0.3 Gyr are clearly inappropri¬ 
ate for these objects as they strongly overproduce the rest- 
frame UV fluxes. A similar result was found by McCarthy 
et al. (2004) and may seem to imply relatively short forma¬ 
tion timescales (e.g. less than 0.3 Gyr). This appears to be 
consistent with the a-enhancement found to be typical of el¬ 
liptical galaxies at z = 0, that also suggests short formation 
timescales (e.g., Thomas, Maraston & Bender 2002; Thomas 
et al. 2005). However, the failure of models with r > 0.3 Gyr 
may not be due to their long formation timescales, but to the 
exponentially declining tail of star-formation. Physical rea¬ 
sons, like the feedback from supernovae explosions or from 
the onset of AGN activity (Granato et al. 2004), suggests that 
in more physically motivated models the SE might be com¬ 
pletely stopped by such feedback. The “truncated" models 
that we considered allow us to explore this possibility, and 
check whether much longer SE timescales are consistent with 
our data. In fact, these “truncated" models produce accept¬ 
able fits to the SEDs even for the case of 2-Gyr lasting star- 
formation. In conclusion, the available portion of the SED 
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B-band 


z-band 


Fig. 10.— Comparison of the surface brightness distribution for object #4950 in the B- and z-band. The same range of fluxes (/x) are plotted for the two 
objects. The peak of the z-band image corresponds to a relative dip in the B-band, that it is instead dominated by two blobs separated by ~ 1 kpc (0.1") from the 
galaxy center. The z-band surface brightness profile clearly shows the central cu™ and extended wings typical of de Vauculeurs galaxies. The faint spiral arms, 
or ring, are not visible in these images. See also Fig.|^and the color inset in Fig.|^for comparison. 


and the presence of the Mguy feature allow us to set strong 
constraints on AgCpass, but not on the previous duration of the 
SF activity. Rest-frame near-IR data from Spitzer could help 
constraining also this latter quantity. 

The hts with r = 0.1 Gyr models allow us to derive lim¬ 
its to the maximum amount of ongoing SFR in these passive 
galaxies, as these models have some residual SF at any time. 
For the 5 objects one gets residual SFR 0. IMQyr^ * with the 
T = 0.1 Gyr models. Even if continuing at this level for a full 
Hubble time to z = 0, this would increase the galaxies masses 
by less than 1%. Of course, we neglect here the possibility of 
completely dust extincted SF that we would not detect anyway 
in the optic al-IR. 

Further constraint on the formation processes of early-type 
galaxies can be obtained from the comparison with theoretical 
simulations. For example Hernandez & Fee (2004) estimated 
the relaxation time-scales of elliptical galaxies formed as the 
result of the merger of two equal size disks, and the Press- 
Schecter probability of such a merger to happen in practice. 
For early-type galaxies with stellar masses similar to the ones 
that we derive, i.e. M* ~ IO^'Mq, they suggest that if such 
objects are found at z 1 . 6 , as we do, it is highly unlikely 
that they have formed through mergers of equal-size spirals. 

6.2. Blue excess early-type galaxies 

We now come back to the two objects with B-band flux 
excess, #4950 and #1025. It is perhaps surprising that Agepass 
measurements for these two objects are not distinguishable 
from the ones of the other redder galaxies. However, unlike 
for the other five galaxies, these two bluest early-type objects 
can only be fitted by the r = 0.3 Gyr models, which seem 
to produce the correct balance between old and younger stars. 
The residual SFRs of the best-fitting t = 0.3 Gyr models is 5- 
10 M 0 yr^' for #4950 and about half of that for #1025. This is 
revealed for #4950 also by the detection of a weak [Oil] A3727 
emission in the K20 survey spectrum. 

These two objects have BzK > —0.2 because of their rela¬ 
tively blue B — z color (see (Fig.0}. Passively evolving galax¬ 
ies shortly after the quenching of SF are expected to show 
similar colors (Daddi et al. 2004b), but on the other hand the 
blue B — z colors may also be due to a secondary burst of star- 
formation, involving a small fraction of the total stellar mass. 


This is most likely the case of the galaxy #4950. Its 
morphology changes strongly from the B-band to the z-band 
(Fig.ia and CoJ. The B-band light is fairly irregular with a 
significant substructure of h-A blobs. It is unlikely that such 
a clumpy structure is a bar, as these are typically more regu¬ 
lar. The blobs visible in the B-band may be instead due to an 
ongoing merger or cannibalism of smaller galaxies or satel¬ 
lites. The z-band light distribution is instead more regular and 
appears typical of a spheroid. A feature looking like spiral 
arms or perhaps a ring is also detected at a distance of about 
5 kpc from the center, that is significantly less luminous than 
the central light both in B and z-bands. From the spectrum 
we And additional evidence for composite stellar populations. 
The best stellar template At is in fact an FO, but best fits with 
A5 and F2 stars are found when limiting the fit to below or 
above 2800A rest frame, respectively. The B-bright clumps 
may have been formed or accreted after the z-band spheroid 
was in place. 

The morphology of object #1025 is instead regular from the 
B-band to the z-band with no obvious changes. Sersic profile 
fits were obtained for this object also in the V- and B-bands, 
finding that it is consistent with an high Sersic index n ~ 4 
in all ACS bands. There is instead a trend of decreasing ef¬ 
fective radius with decreasing wavelength that appears signifi¬ 
cant, possibly due to the presence of a color gradient (Sect.|^. 
In the HST grism spectrum we And that no change is detected 
when fltting templates of stars below or above 2800A rest- 
frame only, and an A3 type star spectrum dominates the light 
consistently in the 2000A< A < 3600A range. There is there¬ 
fore no strong evidence in this object for the coexistence of 
two clearly distinct stellar populations, although the fact that 
a simple SSP model cannot reproduce its SED implies at least 
a somewhat extended SF timescale. This object appears to 
be an early-type galaxy observed very closely to some sig- 
niflcant star-formation event, when early A-type stars (rather 
than F-type stars found in all the others passive objects) were 
still dominating its UV SED. Broadhurst & Bouwens (2000) 
discussed the observational lack of known early-type galaxies 
with A-type dominated stellar populations, our object being 
perhaps the flrst known case. The relatively short lifetime of 
A-type stars, of order of a few hundreds Myr at most, may be 
the actual reason why these objects are rare. 
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The recent star-formation event may actually be the one 
when most of the stars in this spheroid were formed. Ob¬ 
ject #1025 would then represent the closest known link be¬ 
tween passively evolving F-star dominated early-type galax¬ 
ies and the as yet not established nature of the main formation 
event. The relaxed morphology at this early stage would im¬ 
ply a fairly short dynamical relaxation timescale. If the star- 
formation event is a secondary one, involving only a fraction 
of the mass, this objects might instead be similar to the E-l-A 
galaxies found at lower redshift (Dressier et al. 1983). This 
single object in the 1.4 < z < 2 redshift range corresponds, 
however, to a volume density that is 2 to 3 orders of mag¬ 
nitude higher than the local density of E-fA galaxies above 
similar luminosities (Quintero et al. 2004). The fraction of 
E-fA galaxies in the present sample (one out of seven) would 
also be higher then the just 1% that is found locally (Quintero 
et al. 2004). 

7. RESULTS: EVIDENCE FOR A DECLINING NUMBER DENSITY OF 
EARLY-TYPE GALAXIES AT z > 1.4 

An important question to address is how the abundance 
of z ^ 1.4 passive early-type galaxies compares to the local 
value, and to the populations of star-forming galaxies at the 
same redshifts. 

As the majority of our sample galaxies, and all of the class 
"A" ones, are located at z < 2, we first consider the 1.39 < 
z < 2.00 redshift range for this comparison. Over the 12.2 
arcmin^ selection area this redshift range includes a volume 
of ~ 20000 Mpc^. Summing up the galaxy stellar masses 
estimated in Sect.|6lwe obtain a total stellar mass in the range 
of 7-16x IO'^Mq, and a stellar mass density in the range of 
3.4-8.1 X IO^MqMpc^^ for an average redshift < z >= 1.7. 
The 6 early-type galaxies with 1.39 < z < 2.00 correspond to 
a number density of 3.4 x 10-"^ Mpc-3. 

7.1. Evolution from z ^ 2 to z = 0 

The masses of these 6 galaxies are all close to or larger 
than 10^ *M©, and we can assume our sample to be reasonably 
complete above such a mass threshold for passive early-type 
galaxies down to = 23. In the local universe, the num¬ 
ber density and the stellar mass density of passively evolv¬ 
ing galaxies with M* > lO'^M© are about 9 x 10^"*^ Mpc^^ 
and 1.3 X IO^MqMpc^^, respectively (Baldry et al. 2004'^). 
Thus our sample of early type galaxies in 1.39 < z < 2.00 
accounts, at face value, for ^ 35% of the number density of 
M* > 10*^M© passive galaxies in the local universe, and for 
a fraction of ^ 25-60% of the relative stellar mass density. 
Notice that at z = 0 about 75% of all M* > 10**M© galax¬ 
ies are classified as passive early-type galaxies (Baldry et al. 
2004). We therefore infer at z = 1.7 a decrease in the number 
density of massive early-type galaxies by at most of a fac¬ 
tor of 3, if limiting to z < 2. The inferred fractions would 
be further divided by a factor of 2 if using all the volume in 
1.39 <z< 2.47. 

A major limit of these calculations is represented by cos¬ 
mic variance, because of the small volume probed by UDE. 
All our passive galaxies are EROs with R~K > 5.3 (Tabled, 
which are strongly clustered populations (Daddi et al. 2000; 
McCarthy et al. 2001). Clustering can produce strong fluc¬ 
tuations in the number density of objects over small volumes. 

We account for the fact that our pure Salpeter IMF gives stellar masses 
at least 30% larger than for the IMF used in Baldry et al. 2004 


which most likely result in an underestimate of the true densi¬ 
ties as averaged over large volumes (Daddi et al. 2000). If one 
assume that these z > 1.4 galaxies have the correlation length 
of EROs, i.e. ro ~ 10 /i^* Mpc (Daddi et al. 2001; 2003), 
one finds that even at the Icr level the true number density 
(and thus the stellar mass density) of passive galaxies could 
be within half and twice that estimated from our sample, or 
between 20% and 80% of the local value. 

It is clear that similar searches on much larger and indepen¬ 
dent areas would be important to measure the number density 
evolution with some reasonable accuracy. The UDE dataset is 
likely to remain unique in its depth and quality for many years 
to come, but our result confirms that reliable estimates on the 
abundance of z > 1.4 passive early-type galaxies might be ob¬ 
tained by using the BzK photometric criteria only (Eig. 0 
Daddi et al. 2004b), even without complete spectroscopic 
follow-up. In fact, in the UDE the sky density of the brightest 
(Kvega < 20) candidate z > 1.4 early-type galaxies satisfying 
conditions z — K > 2.5 and BzK < —0.2 is not much smaller 
than that found in a 30 times wider field survey (Kong et 
al. in preparation). So, the densities derived over the UDE 
field should not dramatically underestimate the true values. 

7.2. Passive objects at even higher redshifts? 

The highest redshift passive object that we recover is at 
z = 2.47, although with "B" class redshift, while all the others 
lie at z < 2. Therefore, it is important to ask is whether the 
paucity of z > 2 passive galaxies is due to the limiting magni¬ 
tude of the sample (i.e. they exist but are fainter) or to passive 
galaxies getting rarer and rarer. Indeed, some of the galaxies 
in our sample are consistent with having started pure passive 
evolution only at z 2. 

We have used a V /Umax test to try to shed light on the 
issue. Eor each of our objects we compute the maximum 
redshifts at which it would still be part of the sample (i.e. 
with Kab < 23, Kron magnitudes, and z — K > 2.5). We 
use the best fitting Bruzual & Chariot models to compute 
K-corrections, and z = 1.39 as the lower redshift boundary. 
The result is that the six objects would be still in the sample 
at least up to z = 3.3 (when #8238 would drop out) and up 
to z = 4.6 (when finally also #3650 would drop out). With 
< V /Umax >= 0.14 ± 0.04 and a maximum V /Umax = 0.35 
for #1446, an uniform V /Umax distribution is rejected at the 
99.7% level, based on a Kolmogorov-Smirnov test. 

There is therefore strong evidence for evolution, especially 
beyond z = 2. This means that if galaxies with masses and 
SEDs (i.e. ages) similar to the ones that we see at < z >= 1.7 
were present at higher redshifts we would have detected them 
out to z 4 (and with z — K> 2.5). These results under¬ 
score a rapid disappearance of truly passive systems at red¬ 
shifts z ^ 2. This does not mean that massive galaxies where 
not present at earlier epochs, but that if existing they were 
most likely still actively star-forming. Notice that at z = 4 
the Universe is already ^1.6 Gyr old and in principle there 
would have been time to produce passively evolving galaxies 
with ages similar to the ones we see at 1.4 < z < 2.5. 

Of course, clustering can bias also this < V/V,„ax > test. Eor 
example, by chance the UDE may lack a cluster-like structure 
with many similar galaxies at z > 2 or 3. On the other hand, 
a similar < VIVmax > test in the sky contiguous K20 sur¬ 
vey region suggests that the K-band luminous BzK-selected 
starbursts galaxies are evolving positively in number over the 
same redshift range (Daddi et al. 2004b). These BzK star- 
bursts are massive (stellar masses ~ 10**M©) and plausible 


14 


E. Daddi et al. 


candidates for being precursors to early-type galaxies. There¬ 
fore it appears that while the number density of passive galax¬ 
ies is dropping with increasing z, that of their immediate likely 
precursors is somewhat increasing, as expected. It is worth 
emphasizing that these two < V lYmax > tests were performed 
basically on the same sky area, and therefore they should be 
affected by clustering in a similar way. 

7.3. The fraction of z >1.4 stellar mass in passive objects 

It is interesting to estimate what fraction of the mass, or of 
the most massive galaxies, is in passive galaxies at the probed 
redshifts 1.4 < z < 2.5. Galaxies with BzK > —0.2 are in 
most cases star-forming objects at 1.4 < z < 2.5 (Daddi et 
al. 2004b), thus spanning the same redshift range of the de¬ 
tected passive galaxies. We can thus simply compare the sam¬ 
ple of BzK > —0.2 galaxies selected to K^b = 23 to the cur¬ 
rent sample of passive objects to the same limiting magnitude. 
We used the (SED-fit based) /T-band magnitude versus stellar 
mass relation calibrated in Daddi et al. (2004b). Eor the pas¬ 
sive galaxies in our sample this relation provides fairly good 
estimates of the stellar masses, fully consistent with those 
listed in Tabled As star-forming galaxies in 1.4 < z < 2.5 
we consider all BzK > —0.2 sources, i.e. 43 objects when 
excluding the 2 passive "blue" objects and any source with 
X-ray detection as likely an AGN. Of these 6 have estimated 
stellar mass M* ^ 10 * W©, producing a total stellar mass of 
order of 5 x 10* *4/©, or up to 10 x 10* ^Mq for those above an 
estimated stellar mass limit of M* > 0.5 x 10**M©*^. These 
figures are similar to the ones derived above for passive ob¬ 
jects in the same redshift and mass ranges. This supports ear¬ 
lier hints that at 1.4 < z < 2.0 roughly ss 50% of the stellar 
mass is in passively evolving galaxies and roughly « 50% 
is in vigorous star-forming galaxies, for objects with stellar 
masses M* 1 O**M 0 (Cimatti et al. 2004; McCarthy et al. 
2004). By summing together the similar stellar mass densi¬ 
ties contributions from passive and star-forming galaxies at 

1.4 < z < 2 having M* 10**M 0 , one gets fairly close to the 
local value from the SDSS (Baldry et al. 2004). However, 
the uncertainties due to small number statistics and cluster¬ 
ing, and in the estimates of stellar masses are both at least a 
factor of ~ 2. The evolution of most massive galaxies is not 
easily detectable up to z = 2 , and therefore it is arguably not 
very strong. 

7.4. Comparison to KCDM theoretical predictions 

McCarthy et al. (2004) and Cimatti et al. (2004) already 
pointed out that the existence of z > 1.4 passive early-type 
galaxies in such relatively high number is at odd with current 
predictions of semianalytic models of galaxy formation and 
evolution, that predict the passive evolution phase to be estab¬ 
lished much later, if ever. This is similar to the reported fail¬ 
ure to accounts for EROs number counts even at z ~ 1 (e.g., 
Daddi et al. 2000; Eirth et al. 2002). These conclusions are 
supported also by the present findings. 

Eor example, among the early semianalytic models, the 
simulations*^ by Kauffmann et al. (1999) predict a number 
density at z = 1.46 of galaxies with [B — V)rest > 0.6 (Ta¬ 
ble 0 and stellar masses M* > 1 O**M 0 that is about a factor 

The stellar masses estimated for the SF galaxy candidates at 1.4 < z < 

2.5 are not based on magnitudes obtained from surface brightness profile fit¬ 
ting, as for the passive galaxies, and may thus still be somewhat underesti¬ 
mated 

http://www.mpa-garching.mpg.de/GIF/ 
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Fig. 11.— We show the predicted stellar masses versus rest frame (B — V) 
colors of galaxies from the ACDM simulations (Kauffmann et al. 1999). The 
connected large symbols are for our objects, where minimum and maximum 
stellar masses are shown. The early-type galaxies at z < 1.8 are plotted over 
the z = 1.46 simulation, while the ones at z > 1-8 over the z = 2.12 sim¬ 
ulation. The volume of each simulation is 200 times larger than the UDF 
between 1.4 < z < 2.5. 


of 10 lower than recovered here in the UDE, and by z ^ 2 
such objects should have virtually disappeared. Nevertheless, 
it is interesting that these models do predict that at least some 
galaxies exist with the colors and masses of the ones we ob¬ 
serve (Eig. n U . and that a hint for a bimodal color distribu¬ 
tion is already in place in the models at these redshifts. On 
the other hand, there appear to be too few massive galaxies 
in these simulations, either passive or star forming (Daddi et 
al. 2004a,b; Eontana et al. 2004). Among the more recent 
generation of models, those by Somerville et al. (2004) are 
quite successful in predicting the number density of massive 
galaxies at high redshift, but fail to produce the bimodal color 
distribution even at low redshift. Similarly, the hydrodynami- 
cal simulations of Nagamine et al. (2005a;b) appear to match 
the statistic of z = 2 massive galaxies but are less successful 
in reproducing the abundance of z ~ 2 old and passive galax¬ 
ies. There is now general agreement that one of the problems 
with the models is their tendency to sustain star formation in 
massive DM halos all the way to low redshifts, and intense 
theoretical efforts are currently under way in the attempt to 
cope with these discrepancies. A strong AGN feedback ap¬ 
pears to be a viable way for switching off star formation in 
massive galaxies at z ^ 2 (e.g., Granato et al. 2004). 

8 . RESULTS: MORPHOLOGY EVOLUTION OF EARLY TYPE 
GALAXIES TO z > 1.4 

The morphological properties of our sample of z > 1.4 pas¬ 
sive early-type galaxies were compared to those of the corre¬ 
sponding local populations. 

Eig .ll2lshows the B-band rest-frame Kormendy (1977) re¬ 
lation for our objects compared to a local sample of early- 
type galaxies in the Coma cluster (Joergensen et al. 1995) 
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Fig. 12.— The B-band rest-frame Kormendy (1977) relation for the 
1.4 < z < 2.5 early-type galaxies (filled circles with error bars), compared 
to Coma early type galaxies taken from Joergensen et al. (1995) converted 
to the B-band (empty squares), and to 0.7 < z < 1.2 early-type galaxies from 
the K20 survey (crosses, from di Serego Alighieri et al. 2005, in preparation). 
Observed values for our objects are derived from the z-band measurements. 
The average effective brightness < fie> were corrected for cosmological sur¬ 
face brightness dimming and bandpass shift (observed z-band to rest B-band), 
but not for the (relatively uncertain) passive evolution dimming to z = 0. 


and to a sample of massive z ^ 1 early-type galaxies from 
the K20 survey (di Serego Alighieri et al. 2005, in prepa¬ 
ration). The average effective surface brightness < /ig > (see 
e.g., Ziegler et al. 1999 for its definition) has been corrected to 
z = 0 accounting for the (1 -t-z)"*^ cosmological surface bright¬ 
ness dimming. A color term was applied to transform our 
z-band measurements to rest-frame B-band, and similarly the 
Joergensen et al. (1995) r-band measurements were corrected 
using B — r = 1.15 (Ziegler et al. 1999). When not account¬ 
ing for the expected dimming of the B-band rest-frame due 
to the aging of the stellar populations, the z > 1.4 objects are 
largely offsetted from the local Kormendy relation (Fig. [13 
The expected dimming to z = 0 can be computed on the ba¬ 
sis of the inferred best-fitting Bruzual & Chariot (2003) mod¬ 
els, although with large uncertainties given to the range of ac¬ 
ceptable best fitting parameters (especially ages and SF histo¬ 
ries). Typically a dimming by a few magnitudes is derived, 
and these large passive evolution corrections, albeit uncer¬ 
tain, would bring the 7 passive z > 1.4 galaxies in reasonable 
agreement with the local Kormendy relation. 

We notice, however, that 4 out of the 7 objects have very 
small effective radii, rg < 1 kpc, and passive evolution would 
bring them in a region of the — pe diagram where very few 
local galaxies are found. Moreover, local galaxies with Tg < 1 
kpc tend to be much less massive than our z > 1.4 objects. 
Similarly, the stellar surface mass density p* of the 4 galaxies, 
defined as p,* = 44*/ {2r^) is more than a factor ^10 higher 
than that of local galaxies with similar stellar masses (Kauff- 
mann et al. 2004). All in all, the z > 1.4 < 1 kpc objects 

cannot be the progenitors of the z = 0, low-mass galaxies with 
similar effective radii. 


Similarly, about half of z ~ 1 massive early-type galaxies 
from the K20 survey, shown in Fig. [21 have sizes of about 
1 kpc that are smaller then their local massive counterparts. 
Small sizes of K20 early-type galaxies were also inferred by 
Cassata et al. (2005) for the higher redshift sample of 4 early- 
type galaxies at z > 1.5, i.e., = 1-3 kpc, measured on ACS 

z-band data as in this work. Waddington et al. (2002) showed 
that their two z 1.5 radio-selected early-type galaxies have 
too small sizes (3 kpc), as compared to lower redshift samples. 
Stanford et al. (2004) present a rest-frame B-band morpho¬ 
logical analysis of early-type galaxies to z = 1.4 in the Hubble 
Deep Field North based on HSTh-NICMOS, showing indeed a 
trend of significantly smaller sizes of luminous galaxies with 
respect to local massive counterparts (see their Fig. 16). An 
analogous result might be noticed in Fig. 4 of Gebhardt et al. 
(2003). 

Therefore, it appears well established observationally that 
some fraction of high redshift passive early-type galaxies ap¬ 
pear significantly smaller then their likely z = 0 descendants. 
While this is similar to the well known trend of size reduc¬ 
tion with redshift for star-forming galaxies (e.g., Ferguson et 
al. 2004), this effect had not been highlighted and/or dis¬ 
cussed previously for early-type galaxies, if not for the 2 
radio-galaxies of Waddington et al. (2002). It is somehow sur¬ 
prising because one would expect that these massive galaxies, 
that are already undergoing passive evolution, should evolve 
smoothly into local massive ellipticals with no major change 
in their properties if not for aging of the stellar populations. 
The availability of the ultradeep UDF imaging might allow us 
to better constrain and understand this issue. 

We tried to investigate the possible reasons for this effect 
on our sample of z > 1.4 early-type galaxies. The small effec¬ 
tive radii we derive with the Sersic fitting might be due to the 
presence of an unresolved nuclear component, e.g., from an 
AGN. Actually, two of the four objects with small re are de¬ 
tected in X-rays (see Section[3- For the two compact objects 
with the highest S/N ratio in the z-band imaging (#1025 and 
#3650) we attempted fitting their surface brightness profiles 
with two component models, Sersic plus a point-like source. 
The resulting best fit is still dominated by the Sersic compo¬ 
nent, but its re is not significantly increased with respect to 
the one component fits, just a lower Sersic index n is derived. 
However, when forcing to be large, a good fit with ~few 
kpc can be obtained, consistently with the local objects, and 
still with large n values. In this case the nuclear source is 
about 2 magnitudes fainter than the spheroid, in agreement 
with the fact that the optical SEDs and spectra appear domi¬ 
nated by stellar light (Sect. |4]and[3- If this is the reason of 
the small radii, it would imply a widespread nuclear activity, 
or relics of nuclear starbursts, in z ~ 1.4-2 early-type galax¬ 
ies. Close inspection of the inner regions of local early-type 
galaxies does in fact reveal the presence of faint nuclear point 
sources in nearly 50% of the cases (Ravindranath et al. 2001). 
Those nuclear sources might have been significantly more lu¬ 
minous/frequent at high redshifts. 

If instead the small values are not (or not always) due to 
the presence of nuclear emission, two possibilities may help to 
solve the puzzle, and reach consistency with the local galaxy 
populations: either a strong morphological K-correction or 
the physical sizes will increase with time. Larger sizes at 
z = 0 would move the galaxies roughly along the Kormendy 
relation by decreasing the average effective surface bright¬ 
ness, and would bring our objects positions in Fig.^jin better 
agreement with local ones. 
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Fig. 13.— The (i — z) color maps: darker regions correspond to bluer i~z 
colors. Only pixels with S/N> 3 in both i and z bands were considered. No 
attempt was made to account for the different PSF between the two bands, so 
that the bluer (i.e. darker in the images) cores of most objects are likely the 
result of PSF differences. Obvious exceptions are object #4950 where 2 blue 
blobs are clearly visible, and #1025 that shows a strong color gradient and 
blue core. 


A possible mechanism for size evolution would be substan¬ 
tial satellite engulfment by the massive red galaxies as they 
evolve to lower redshifts, which would imply an additional 
growth of their stellar mass. Object #4950 might be an exam¬ 
ple of this process in act. Confirmation of this scenario would 
require the detection of large concentrations of satellite galax¬ 
ies in the vicinity of the passive z > 1.4 galaxies. This option 
might be tested in the future by using the ultradeep UDF data. 

In relation to the morphological K-correction option, we 
note that blue cores are not uncommon already among 0.5 < 
z < 1 ellipticals (Menanteau et al. 2001, 2004), and that the 
observed z-band corresponds to 2600-3800A in the rest frame 
for 1.4 < z < 2.5 (so somewhat bluer than the 4300A of the B- 
band). The most direct way to check this option would be to 
obtain Sersic profile fits in the NICMOS FI lOW and F160W 
bands, bracketing the rest-frame Z?-band for 1.4 < z < 2.5. 
However, this is technically difficult, mainly due to the large 
0.2" pixel size of NICMOS (a factor of 4 larger than the ACS 
pixels) corresponding to ^ 2 kpc at z ~ 2. Attempting such 
fits is beyond the scope of the present work. What we checked 
instead was the trend of with wavelength, within the range 
covered by the ACS data. We constructed i — z color maps 
of all our z > 1.4 galaxies (Fig. [nil, but among the 4 objects 
with re < 1 kpc only #1025 shows a quite prominent blue core. 
This objects also exhibits a strong trend of rg with wavelength, 
as shown in Fig. d while a Sersic profile consistent with 
n ^ 4 is found in all the ACS bands. An extrapolation would 
give re larger by 50% in the rest-frame B-band with respect 
to that measured in the observed z-band. Still, this increase 
would not suffice to bring this object at z = 0 on the pop¬ 
ulated part of the Kormendy relation, as that would need an 



Rest-Irame wavelength (A) 


Fig. 14.— The dependence of the effective radius on wavelength for object 
#1025 at z = 1.73. The observed trend would imply about 50% larger at 
rest-frame B-band (4350A) with respect to observed z-band (or rest-frame 
3350A), but still ~ 1 kpc only. 


increase by a factor ^ 3. Yet, we cannot presently rule-out the 
possibility that rg may change more steeply with wavelength, 
especially beyond the 4000A break. We finally note that small 
re values in the rest-frame B band are common among early- 
type galaxies to z = 1.4 (Stanford et al. 2004), which outlines 
the existence of this problem also when using B-band rest- 
frame measurements. A final assessment of why some high-z 
ellipticals have such small effective radii and how they relate 
to local E/SO galaxies is left to future investigations. 

9 . RESULTS: THE X-RAY PROPERTIES OF z > 1.4 EARLY-TYPE 
GALAXIES 

Two of the 7 z > 1.4 early-type galaxies discussed here are 
detected in the 1 Msec Chandra X-ray observation (Giacconi 
et al. 2002), and their properties are reported in Table|4] The 
two sources are both very hard (NB > 0.4) and both of them 
were discussed also by Padovani et al. (2004) as possible 
QS02 candidates*^. 

Object #1446 is close to the detection limit with ~ 20 net 
counts, almost all from the hard band (HR> 0.42), and we 
estimate its luminosity to be F 2 -iokeV = 2.2 x 10 ^^ erg s^*. 
Object #1025 is detected with ^ 90 net counts, thus allowing 
a rough spectral analysis. Fig. ini shows the X-ray spectrum 
of the object along with an absorbed power-law at the source 
redshift. The best fit rest frame column density id Ah=5 x 10^3 
cm^^; and the corresponding unabsorbed X-ray luminosity is 
L 2 -iokeV= 8 x 10^^ erg s^*. Given the high column density 
and the X-ray luminosity close to 10^^ erg s^* this object 
can be classified as an X-ray QS02. It is intriguing that the 
optical spectrum appears to be dominated by stellar light, that 
may imply that the obscuring torus is blocking also most of 

The luminosities derived by Padovani et al. (2004) are larger than those 
inferred here given that these authors put these objects at slightly higher red- 
shifts. 
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Table 4. X-ray emission properties 


ID 

ID(G02) 

ID(A03) 

z 

^soft 

Fhard 

Nh 

G-lOteV 

X/O 

X/K 





cgs 

Cgs 


Cgs 



1025 

256 

255 

1.73 

0.9 X 10-‘® 

4.1 X 10^*^ 

~ 5 X 10^3 

8 X 10“*^ 

12 

2.5 

1446 

605 

243 

2.47 

<0.4x lO^'*’ 

4.6 X 10-‘® 

(HR> 0.42) 

2 X 10“*^ 

5 

1 


the emission in the optical and near-IR domains. An excess 
emission with respect to the power-law continuum is present 
at ~ 2.3 — 2.4 keV. This would be consistent with the presence 
of a redshifted 6.4 keV Fe Ka emission line. Indeed, when a 
Gaussian line is added to the model spectrum, the fit improves 
and suggests a redshift z = 1.71, in excellent agreement with 
the redshift measured in the optical/near-IR, and a rest-frame 
line equivalent width of ~ 350 eV, similar to what found in 
obscured X-ray sources (e.g., Brusa, Gilli & Comastri 2005b; 
Maccacaro et al. 2005). 

Both sources have large X-ray to optical flux ratios (X/O). 
Coupled with the analysis of Mignoli et al. (2004), that sug¬ 
gest that most of the hosts of hard X-ray AGN with large X/O 
and R-K> 5 are spheroids at z > 1, this imply a close connec¬ 
tion between hard X-ray sources and high-redshift early-type 
galaxies. 

Although based on two objects only, the fraction of AGNs 
in z ~ 2 these early-type galaxies is at face value ^ 30%. This 
is larger than that observed for the early-type sample of z ~ 1 
EROs discussed by Brusa et al. (2002; no X-ray emission was 
revealed in a sample of 8 sources), and much larger than the 
~ 1% value observed locally. This is similar to the fraction of 
AGN among EROs at the probed K magnitudes (Alexander et 
al. 2002; Brusa et al. 2005a) underlinyng again the possible 
close connection between the formation of AGN and early- 
type galaxies. 



Fig. 15.—The Xray spectrum of #1025 from Chandra, fitted with a F = 2.0 
power law and absorbing column density of Nh = 5 x 10^^ cm“^. The excess 
flux at 2.4 keV is consistent with the 6.4 keV Fe line, redshifted at z = 1.73. 
The presence of such a strong Fe emission is reasonable (Brusa et al. 2005b), 
given also the very hard X-ray spectrum and high absorption column density 
of the source. 


10 . SUMMARY AND CONCLUSIONS 

Based on the BzK two-color selection criterion meant to 
identify z > 1.4 passively evolving galaxies, a sample of seven 
objects has been identified in the UDE for which ACS grism 
spectra confirm both the high redshift (1.4 z 2.5) and the 
current absence (or very low) ongoing star formation. The 
analysis conducted over the ultradeep ACS imaging further 
shows that the objects are morphologically early-type galax¬ 
ies. 

The redshifts are derived from the identification of the spec¬ 
tral feature at the rest-frame 2640 < A < 2850 A due mainly 
to Mgl and Mgll absorptions in late A- and E-type stars, and 
which shows up prominently only in synthetic stellar popu¬ 
lations with no (or very low) star formation in the last ^0.5 
Gyr. 

The SEDs and ACS grism spectra of 5 of the objects are 
best reproduced by a star formation history initiated ^1.5 
Gyr before, and completely discontinued over the last ^ 0.5 — 
1.5 Gyr, implying that the passive evolution started at z — 2 — 
3.5. In the other cases a similar age constraint is inferred, but 
the spectra are best reproduced by an exponentially declining 
SER, with T ~ 0.3 Gyr. 

The overall spectral energy distributions of the objects in¬ 
dicate stellar masses in excess of ^ 10^ *M© and we infer 
that the comoving volume density of such massive, early-type 
galaxies at < z >— 1.7 is about 1 /3 of the local value. The un¬ 
certainty in the derived fraction is dominated by cosmic vari¬ 
ance, given the small field covered by the UDE. Allowing for 
the expected clustering of these galaxies we infer that at the 
1(7 level such fraction could be as low as ~ 20% or as high 
as ~ 80%. Clearly, the exploration of wider fields is neces¬ 
sary to accurately pinpoint the actual decrease with redshifts 
of the number density of massive galaxies which are passively 
evolving. However, a V /Umax test indicates that beyond z ~ 2 
this number density is likely to drop very rapidly. 

One intriguing aspect of the present findings is the quite 
small effective radii derived from the ACS z-band images, 
which in 4 out of 7 cases are less than 1 kpc. Given their 
high mass, the passive evolution of such very compact ob¬ 
jects will bring them in a region of the Kormendy relation (or, 
equivalently, of the fundamental plane) which is depopulated 
at z = 0. We discuss various possibilities to explain this ap¬ 
parent paradox, including a morphological K-correction (e.g., 
due to the possible presence of blue cores in these galaxies), or 
of an AGN point-like source biasing the r,, measurements, or 
eventually some evolutionary effects such as if the observed 
galaxies were still subject to further growth at lower redshifts. 
While we mention hints favoring one or another option, the 
existing data do not allow to reach any firm conclusion on 
this specific issue. 

Einally, we point out that while in the local universe most 
of the most massive galaxies are passively evolving giant el¬ 
lipticals, the present data - combined with previous evidence 
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in a partly overlapping field (Daddi et al. 2004a;b) - indi¬ 
cates that by z ^ 2 passively-evolving and seemingly vigorous 
starbursts galaxies occur in comparable numbers among most 
massive galaxies. This remains a major challenge for the cur¬ 
rent theoretical simulations of galaxy formation to reproduce. 
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